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In this work, the electronic properties of phosphorene nanoribbons with different width and edge 
configurations are studied by using density functional theory. It is found that the armchair phosphorene 
nanoribbons are semiconducting while the zigzag nanoribbons are metallic. The band gaps of armchair 
nanoribbons decrease monotonically with increasing ribbon width. By passivating the edge phosphorus 
atoms with hydrogen, the zigzag series also become semiconducting, while the armchair series exhibit a 
larger band gap than their pristine counterpart. The electronic transport properties of these phosphorene 
nanoribbons are then investigated using Boltzmann theory and relaxation time approximation. We find 
that all the semiconducting nanoribbons exhibit very large values of Seebeck coefficient and can be further 
enhanced by hydrogen passivation at the edge. Taking pristine armchair nanoribbons and 
hydrogen-passivated zigzag naoribbons with width N=7, 8, 9 as examples, we calculate the lattice thermal 
conductivity with the help of phonon Boltzmann transport equation and evaluate the width-dependent 
thermoelectric performance. Due to significantly enhanced Seebeck coefficient and decreased thermal 
conductivity, we find that at least one type of phosphorene nanoribbons can be optimized to exhibit very 
high figure of merit (ZT values) at room temperature, which suggests their appealing thermoelectric 
applications. 



Recently, the black phosphorous has become the focus of science community due to the synthesis of its 
two-dimensional form, namely, the phosphorene. Preliminary but exciting results indicate that phosphor- 
ene has potential application in nanoelectronics and optoelectronics^'^'^. As known, black phosphorus (BP) 
is the most stable phosphorus allotrope under normal condition and owns a direct band gap of about 0.3 eV^'^, 
which can be easily turned into bothp-type and n-type configuration. What's more, BP was found to have high 
carrier mobility in the order of 10^ cm^/V/s at low temperature^'^. Such superior electronic and transport prop- 
erties suggest that BP may have potential thermoelectric applications. Indeed, recent theoretical work of Lv et al 
showed that the power factor of BP can reach as high as 1 18.4 |iWcm"^K"^ by appropriate doping. However, due 
to relatively higher lattice thermal conductivity (12.1 W/mK), the maximum ZT value of BP can only reach 0.22 at 
room temperature At the same time, Qin et al investigated the strain effect on the thermoelectric properties of 
BP^\ They found that applying strain is a practical way to enhance the thermoelectric performance of BP, and the 
largest ZT value of 0.87 can be achieved at 800 K with a tensile strain of 7%. Although these works suggest the 
possibility to use BP as a new kind of thermoelectric material, the ZT values obtained are very small and can not 
competitive with typical thermoelectric materials such as Bi2Te3 and CoSbs. Another efficient method to improve 
the thermoelectric performance is using low- dimensional structures. As theoretically predicated by Hicks and 
Dresselhaus^^'^^, low- dimensional systems could have dramatically larger ZT values than the corresponding bulk 
materials because of decreased thermal conductivity caused by phonon boundary scattering and improved power 
factors on account of quantum confinement. Indeed, Fei et al. showed that the ZT value of two-dimensional 
phosphorene (similar to graphene) is greater than 1.0 at room temperature and can reach 2.5 at 500 K^'*, which 
can be traced to a large ration of electrical conductance to thermal conductance with orthogonally preferred 
conducting directions. Lv et al. demonstrated that strain induced band convergence in two-dimensional phos- 
phorene can significantly increase the Seebeck coefficient and electrical conductivity at the same time, which leads 
to a room temperature ZT value of about 2.1 in the armchair direction^^. If the dimension can be further reduced, 
say, considering the so-called phosphorene nanoribbons (in analogy with graphene nanoribbons), one may 
expect even better thermoelectric performance. It is thus of both fundamental and practical interests to discuss 
the thermoelectric properties of the one- dimensional phosphorene nanoribbons. 
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In this work, we show by first-principles calculations that phos- 
phorene nanoribbons (PNRs) with zigzag (ZPNRs) and armchair 
(APNRs) edges have distinct electronic properties, which can be 
further tuned by hydrogen passivation. Using the Boltzmann theory 
for both electrons and phonons, the transport coefficients of the 
PNRs are calculated. As a quick understanding of the thermoelectric 
performance of such kind of low- dimensional system, the ZT value of 
APNRs and H-passivated ZPNRs with several typical width are eval- 
uated. Our results indicate that the ZT values of APNRs can be 
optimized to as high as 6.4 at room temperature, which suggests that 
phosphorene nanoribbons could be very promising candidates for 
high-performance thermoelectric applications. 

The PNRs can be obtained by cutting a monolayer phosphorene 
along armchair or zigzag directions. Following the conventional 
notation for graphene nanoribbon^^, the armchair phosphorene 
nanoribbons (APNRs) or zigzag phosphorene nanoribbons 
(ZPNRs) can be identified by the number of dimer lines or the zigzag 
chains across the ribbon width and are labeled as AT- APNRs or N- 
ZPNRs, respectively. Figure 1(a) and 1(c) show the ball -and- stick 
model of AT- APNRs and AT-ZPNRs, respectively. Here we consider 
N=7 ~ 12 for both armchair and zigzag nanoribbons, and the width 
varies from about 10 A to 26 A. Upon structure relaxations, we see 
from Figure 1(b) and 1(d) there are some edge reconstructions for 
both APNRs and ZPNRs. Take N=9 as an example. For the APNR, 
we find that the bond length between the edge atoms on the trans- 
lation direction decreases from 2.26 A to 2.06 A. The corresponding 
bond angles at edge increase from 96.9° to 111.0° for a, and from 
103.7° to 119.5° for p. In the case of ZPNR, the bond length at the 
edge decreases from 2.23 A to 2.14 A, and the corresponding edge 
angles y and 6 increase from 96.9° to 100.7°, and from 103.7° to 
108.9°, respectively. Similar edge reconstructions have been also 
found in previous calculations^^. It should be mentioned that the 
reconstructed motifs at the opposite edges are symmetrically 
arranged, which is different from those found in the BiSb nanorib- 
bon^^ with similar atomic configurations. 

To investigate the relative stability of pristine phosphorene nanor- 
ibbons, we have computed the formation energy of both APNRs and 
ZPNRs. Here the formation energy EpE is defined as: 



EpE = (Ep 



-nEp)/n, 



(1) 



where Ep^ji is the calculated total energy of a given PNR, Ep is the 
energy per atom of a phosphorene sheet, and n is the number of P 
atoms in the PNR. Figure 2(a) plots the formation energies of APNRs 




(b) (d) 

Figure 1 | The ball- and- stick model of armchair phosphorene 
nanoribbon: (a) is top-view, and (b) is side-view. The zigzag counterpart is 
plotted in (c) and (d), respectively. The unit cell is outlined by dash lines 
and the arrow indicates the translation direction of nanoribbons. 



and ZPNRs as a function of the ribbon width. The small but positive 
formation energies indicate that the one-dimensional phosphorene 
nanoribbons are energetically less favorable than the two- 
dimensional phosphorene. Compared with the well-known gra- 
phene nanoribbons^^, we find that both armchair and zigzag phos- 
phorene nanoribbons have much smaller formation energies, 
suggesting that it is quite accessible to synthesize the PNRs experi- 
mentally. On the other hand, we see that the ZPNRs are energetically 
more favorable than the APNRs. The formation energies of both 
systems decrease monotonically with increasing ribbons width and 
converge to zero, indicating that the wider nanoribbon has relatively 
higher stability. To further confirm the stability of the phosphorene 
nanoribbons, we have performed ab-initio molecular dynamics 
(MD). The starting temperature is 300 K and a microcanonical 
ensemble has been used. The MD runs for 2000 steps with a time 
step of 0.5 fs. Taking 7-APNR and 7-ZPNR as an example, we show 
in Figure 2(b) the nearest P-P distance at the edge in the whole MD 
running time. We see there are only slight fluctuations around the P- 
P bond length of 2.06 A for 7-APNR and 2.14 A for 7-ZPNR, which 
indicates that the phosphorene nanoribbons considered in our work 
are rather stable. 

It is known that monolayer phosphorene is semiconducting with a 
direct band gap of 1.51 eV^ By cutting the monolayer into nanor- 
ibbons, however, we find that the ZPNRs become metallic while 
APNRs remain semiconducting but exhibit an indirect band gap. 
As indicated in Figure 3, the band gaps of APNRs decrease with 
increasing ribbon width, which can be attributed to the well-known 
quantum confinement effect^°. If the edge phosphorus atoms are 
passivated by hydrogen atoms, our calculated results indicate that 
the above-mentioned edge reconstructions disappear, and the band 
gaps of H-saturated APNRs become direct. Moreover, there are obvi- 
ous increases of the corresponding band gaps (Fig. 3). Such effect is 
more pronounced for the ZPNRs, where a transition from metal to 
indirect band gap semiconductor is observed, and the H-passivated 
ZPNRs exhibit the largest band gaps among the four kinds of phos- 
phorene nanoribbons. Our results are consistent with previous 
works^^'^^ and further confirm the reliability of our calculations. 

We next consider the electronic transport properties of these 
nanoribbons by using the semiclassical Boltzmann theory and 
rigid-band approach^^. Within this method, the chemical potential 
pi indicates the doping level (or carrier concentration) of the system, 
and the positive and negative pi correspond to n-type and p-type, 
respectively. The optimal carrier concentration can be obtained by 
integrating the density of states (DOS) of the system from the desired 
chemical potential to the Fermi level (^^=0). The Seebeck coefficient 
S and the electrical conductivity a are given by integration of the so- 
called transport distribution H(£) = v-t-: 

' ^ k k k 



-^1 



de 



E{s) 



ksT 



^j..(-|)H,.,, 



(2) 



(3) 



where /o is the equilibrium Fermi function, ks is the Boltzmann's 

constant, and p is the chemical potential, and is the group 

k k 

velocity and relaxation time at state k, respectively. 

Figure 4(a) and 4(b) show the calculated Seebeck coefficients S for 
the APNRs and H-passivated APNRs as a function of chemical 
potential at room temperature, respectively. We see that by optim- 
izing the carrier concentration, the phosphorene nanoribbons can 
exhibit very large value of Seebeck coefficients. As the definition of 
cross-sectional area has some arbitrariness for low-dimensional sys- 
tems such as our phosphorene nanoribbons, in the following we will 
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Figure 2 | (a) The calculated formation energy of APNRs and ZPNRs as a function of the ribbon width, and (b) the ab-initio MD results of the nearest P-P 
distance at the edges for the 7-APNR and 7-ZPNR at room temperature. 



calculate the carrier concentration with respect to the interlayer dis- 
tance of bulk phosphorus (3.07 A) and the real ribbon width. For 
example, the pristine 7-APNR has a maximum Seebeck coefficient of 
1.36 mV/K at ^=0.056 eV. Upon edge passivation, the maximum 
Seebeck coefficient can further enhanced 1.75 mV/K, which is sig- 
nificantly higher than most of previous results^^"^^. As all the inves- 
tigated ZPNRs in their pristine form are metallic, we see from 
Figure 4(c) that their Seebeck coefficients are very small around 
the Fermi level. However, we observe a much higher Seebeck coef- 
ficient for the H-passivated ZPNRs (Figure 4(d)), which can be opti- 
mized to 2.6 mV/K for the ZPNRs with width of N=7. Such giant 
Seebeck coefficients of phosphorene nanoribbons are very beneficial 
for their thermoelectric applications, and we will come back to this 
point later. On the other hand, we find from Figure 4 that the Seebeck 
coefficients of three kinds of semiconducting nanoribbons exhibit 
obvious edge dependence, which is the largest for the H-passivated 
ZPNRs, the smallest for the pristine APNRs, with H-passivated 
APNRs in between. Such an order of Seebeck coefficients is a con- 
sequence of their distinct energy gaps. As shown in Figure 5, the 
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Figure 3 | The calculated energy band gaps for phosphorene nanoribbons 
as a function of width. 



Seebeck coefficient at the peak shows linear dependence with the 
band gap and has no trend towards to saturation. This implies that 
the structures which have larger band gaps would have higher 
Seebeck coefficients. Similar relation has also been found in the 
graphene/h-BN superlattice nanoribbons^^. 

Within the semiclassical Boltzmann theory, the electrical conduc- 
tivity (7 can only be calculated with the relaxation time t inserted as a 
parameter. This means that what we actually obtained is ct/t. To 
figure out the electronic relaxation time of phosphorene nanorib- 
bons, we apply the deformation potential (DP) theory proposed by 
Bardeen and Shockley^^, where the relaxation time t of ID system at 
temperature T can be expressed as: 



{2nkBTf^\m*\'^^Ei^' 



(4) 



In this formula, and are the effective mass and carrier mobility 
along the transport direction, respectively. The deformation poten- 
tial constant Ei = A£/(A///) is determined by changing the lattice 
constant / along the transport direction. The elastic modular can be 
obtained by C = [d^E/dd^]/ly where E is the total energy of the system, 
3 is the applied uniaxial strain along the ribbon direction. 

As the H-passivated zigzag ribbon has the largest Seebeck coef- 
ficient while it is the smallest for the pristine armchair ribbon, for 
simplicity, in the following we take the APNRs and H-passivated 
ZPNRs with width N=7, 8, 9 as examples to discuss the thermoelec- 
tric performance of phosphorene nanoribbons. Table I lists the cal- 
culated room temperature carrier mobility fic^ the relaxation time (t), 
and all the related parameters indicated in Eq. (4) for the pristine 
armchair and H-passivated zigzag ribbons. Note that the calculated 
mobility for the electron can reach as high as 740 cm^/V/s at room 
temperature, which is higher than the experiment values of two- 
dimensional phosphorene FET^'^ and monolayer MoS2^^. The relaxa- 
tion time for electron is higher than that for hole in APNRs, but it is 
just reversed for the H-passivated ZPNRs. 

Figure 6 (a) shows the calculated electrical conductivity a of 7, 8, 9- 
APNR as a function of chemical potential fi at 300 K. We see that cr of 
n-type is usually larger than that ofp-type, especially for the 9-APNR 
which has a much larger electron relaxation time. For both 7-APNR 
and 8-APNR, we see that cr increases gradually with increasing abso- 
lute fi. However, the electrical conductivity is very small at the chem- 
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Figure 4 | The calculated Seebeck coefficient of (a) APNRs, (b) H-saturated APNRs, (c) ZPNRs, and (d) H-saturated ZPNRs as a function of chemical 
potential at 300 K. 



ical potential where the Seebeck coefficient approach its peak value 
(See Fig. 4(a)). This implies that there should be a certain comprom- 
ise between the Seebeck coefficient and the electrical conductivity. 
Indeed, we see from the Figure 6(b) that the maximum power factor 
(denoted by S^g) appears at a particular chemical potential where 
neither the Seebeck coefficient nor the electrical conductivity reach 
the maximum. Among the investigated 7, 8, 9-APNR, we find that 
the 9-APNR exhibit the smallest Seebeck coefficient but the largest 
power factor (n-type), which results from its much larger electrical 
conductivity as discussed above. The calculated electronic thermal 
conductivity (k;^) shown in Figure 6(c) has the same behavior as that 
of electrical conductivity since is calculated by the Wiedemann- 
Franz law^° 
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where the Lorentz number is set to 2.45 X 10"^ WQR-^^^-^l 

To estimate the lattice thermal conductivity {Kp) of phosphorene 
nanoribbons, the phonon Boltzmann transport equation for the pho- 
nons are adopted^^"^^. The calculated lattice thermal conductivity of 
7, 8, 9-APNR and H-passivated 7, 8, 9-ZPNR are Ust in the Table II. 
We find that the lattice thermal conductivity of both series increase 
with the ribbon width, which is consistent with those found in gra- 
phene nanoribbons and silicene nanoribbons^^'^^. In the case of 
APNRs, the lattice thermal conductivity are much smaller than that 
of bulk value (12.1 W/mK)^^, implying that pristine armchair phos- 
phorene nanoribbons may have better thermoelectric performance 
than their bulk counterpart. The reduced lattice thermal conductivity 
can be attributed to the fact that phonons in the armchair nanorib- 
bons are different from those in bulk phosphorus because the dis- 
persion relations will be significantly modified due to confinements 
in two directions. The phonon-phonon interactions can be changed 
and boundary scattering is much stronger in nanoribbons than in 
bulk system. Compared with the electronic lattice thermal conduc- 
tivity (see the inset of Fig. 6(c)), the lattice part of APNRs is relatively 
higher in a certain chemical potential range around the Fermi level. 
However, the electronic part becomes dominative and increases 
rapidly beyond this range. 

With all the transport coefficients available to us, we can now 
evaluate the ZT value of phosphorene nanoribbons, which is given 
by: 
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Figure 5 | The room temperature Seebeck coefficients (peak value) as a 
function of band gaps for APNRs, H-saturated APNRs, and H-saturated 
ZPNRs. 



Figure 6(d) shows the calculated ZT value of 7, 8, 9-APNR as a 
function of chemical potential at room temperature. We see that 
n-type APNRs exhibit larger ZT value than the p -type systems, which 
is mainly caused by relatively larger electrical conductivity and 
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Table 1 The deformation potential E-\, stretching modulus C, effective mass m*, mobility fi^, 
along periodic direction in the 7, 8, 9-APNR and H-passivated 7, 8, 9-ZPNR at 300 K 

N-PNRs Carriers Ei (eV) C(eV/A) m* (mj 


and relaxation time t 

fi, (cm^V-^s-i) 


for hole and electron 

T(fs) 


7 APMP 
/ -Mr INK 


hole 


— 1 /I9n 


/ .00/ 


n AA 


1 no 0 

1 KJZ.Z 


zO.o 




electron 


-0.826 


7.667 


0.65 


1 00.0 


AO 0 

Oz.y 


8-APNR 


hole 


-2.106 


9.776 


0.29 


1 1 1 T 
III./ 


1 o.D 




electron 


-1.087 


9.776 


0.63 


1 Jz.4 


/(T 0 
4/ .Z 


9-APNR 


hole 


-2.569 


12.033 


0.23 


1 oU./ 


1 T 0 
1 /.Z 




electron 


-0.512 


12.033 


0.62 


/4U.U 


OAO ^ 


7-H-ZPNR 


hole 


-0.666 


105.179 


3.67 


268.7 


560.3 




electron 


3.372 


105.179 


1.29 


50.1 


36.8 


8-H-ZPNR 


hole 


-0.801 


120.705 


3.69 


211.4 


443.4 




electron 


3.492 


120.705 


1.30 


52.9 


39.2 


9-H-ZPNR 


hole 


-0.896 


133.838 


3.60 


194.3 


397.5 




electron 


3.351 


133.838 


1.30 


64.3 


47.3 



Seebeck coefficient of the former. Moreover, the peak ZT values ofp- 
type systems decrease with increasing ribbon width, but the behavior 
is reversed for the n-type systems. Among the investigated APNRs, 
the 7-APNR shows the largest p-type ZT value of 2.3 at chemical 
potential fi= —0.36 eV, while the 9-APNR exhibits the largest n-type 
ZT value of 6.4 with ili= 0.21 eV. Such ZT values significantly 
exceeds the performance of most laboratory results reported and 
are comparable to the efficiency of traditional energy conversion 
method, which makes armchair phosphorene nanoribbon a very 
promising candidate for thermoelectric applications. As summarized 
in Table II, the effectiveness of the phosphorene nanoribbon is 
mainly due to its very large Seebeck coefficient and much lower 
thermal conductivity (especially the lattice part) compared with 
the bulk phosphorous 



We next consider the H-passivated ZPNRs which have the largest 
Seebeck coefficients among the four nanoribbon series (see Fig. 4). 
The calculated electronic conductivity cr of 7, 8, 9-ZPNR as a function 
of chemical potential fi is shown in Figure 7(a). Compared with the 
APNRs, the obvious difference is that the p-type a of H-passivated 
ZPNRs is obviously larger than that of n-type, which also leads to 
much larger power factor S^g of the former (see Figure 7(b)). 
However, we see from Table II that H-passivated ZPNRs have 
obviously larger lattice thermal conductivity compared with 
APNRs. Similar to those observed in graphene nanoribbons^^ and 
silicene nanoribbons^^, such strong anisotropy of thermal conduc- 
tivity is mainly caused by different phonon scattering rates at the 
zigzag and armchair edges^^, and low frequency phonon bands of 
zigzag nanoribbons are in general more dispersive than those of the 




0.0 

//(eV) 

Figure 6 | The calculated electronic transport coefficients at 300 K as a function of chemical potential for the 7, 8, 9-APNR: (a) the electrical 
conductivity, (b) the power factor, (c) the electronic thermal conductivity, and (d) the ZT value. The straight lines in the inset of (c) indicate the lattice 
thermal conductivity of 7, 8, 9-APNR. 
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Table II Maximum ZT values and corresponding electron transport coefficients at optimized chemical potential (or carrier concentration 
n) for the 7, 8, 9-APNR and H-passivated 7, 8, 9-ZPNR at room temperature. Note for the one-dimensional PNR, the carrier concentration 
is calculated with respect to the interlayer distance of the bulk phosphorus (3.07 A) and the real ribbon width 

N^PNRs Carrier type (eV) n (cm"^) 5 (^iV/K) a(]0'S/m) S^cr (W/mK^) Ke (W/mK) K;p(W/mK) IT 


7-APNR 


p-type 


-0.36 


5.22 X 10^9 


343.1 


0.86 


1.01 X 10-2 


0.63 


0.67 


2.3 




n-type 


0.35 


3.71 X 10^9 


-431.3 


1.34 


2.49 X 10-2 


0.98 


34.64 


4.5 


8-APNR 


f>type 


-0.28 


4.29 X 10^9 


316.6 


0.79 


7.93 X 10-2 


0.58 


0.75 


1.8 




n-type 


0.27 


4.77 X 10^9 


-435.7 


1.59 


3.02 X 10-2 


1.17 


0.75 


4.7 


9-APNR 


f>type 


-0.25 


5.67X 10^9 


270.8 


1.20 


8.80 X 10-2 


0.88 


1.67 


1.0 




n-type 


0.21 


2.32 X 10^9 


-492.4 


4.13 


1.00 X 10-^ 


3.03 


1.67 


6.4 


7-H-ZPNR 


p-type 


-0.79 


5.12 X 102° 


244.2 


17.4 


1.04 X 10-^ 


12.78 


33.88 


0.7 




n-type 


0.81 


2.90 X 102° 


-193.2 


2.16 


8.07 X 10-2 


1.60 


33.88 


0.1 


8-H-ZPNR 


p-type 


-0.74 


4.61 X 102° 


231.5 


13.5 


7.23 X 10-2 


9.91 


34.64 


0.5 




n-type 


0.76 


2.34 X 102° 


-193.5 


1.97 


7.38 X 10-2 


1.45 


34.74 


0.1 


9-H-ZPNR 


f>type 


-0.69 


4.75 X 102° 


222.4 


12.4 


6.1 1 X 10-2 


9.08 


37.09 


0.4 




n-type 


0.70 


2.24 X 102° 


-193.0 


2.16 


8.04 X 10-2 


1.59 


37.09 


0.1 



armchair nanoribbons^^'^^. As a result, the optimized ZT values of H- 
passivated ZPNRs are rather small, which is only 0.7 for the case ofp- 
type and 0.1 for the n-type system. 

In summary, we have demonstrated that at least one type of phos- 
phorene nanoribbons (APNRs) could be optimized to exhibit very 
good thermoelectric performance with a ZT value as high as 6.4 at 
room temperature. Our theoretical calculations are self- consistent 
within the framework of density functional theory and Boltzmann 
theory, and there are no adjustable parameters. We want to mention 
that the phosphorous thin films and few-layer or even monolayer 
phosphorene were recently fabricated by using mechanical exfo- 
liation method^"^. It is reasonable to expect that by cutting monolayer 
phosphorene or by patterning epitaxially grown phosphorene, one 
can obtain the phosphorene nanoribbons with controllable width 
and/or edge, and thus realize their possible thermoelectric applications. 



Methods 

The electronic structure properties of PNRs are performed by using the first- 
principles plane-wave pseudopotential formulation"^""*^ as implemented in the 
Vienna ab-initio Simulation Package (VASP)*^. The exchange-correlation energy is in 
the form of Perdew-Burke-Ernzerhof (PBE)"^* with generalized gradient approxi- 
mation (GGA). The cutoff energy for the wave function is set to 400 eV. For the 
geometry optimization, a 1 XIX 15 Monkhorse-Packic-meshes'^^ is adopted for the 
Brillouin zone integration. The atomic positions are fully relaxed until the force on all 
atoms become less than 0.05 eV/A. A vacuum distance of 14 A is used for the 
directions of both width and thicknesses, so that the nanoribbons can be treated as 
independent entities. Based on the calculated energy band structure, the electronic 
transport properties are determined by using the semiclassical Boltzmann theory 
within the relaxation time approximation"^*^. Such a method has been successfully used 
to predict some known thermoelectric materials, and the theoretical calculations 
agree well with the experimental results^^"^*^. In order to obtain reliable transport 
coefficients, a denser Monkhorse-Pack ic-mesh up to 125 points in the irreducible 
Brillouin zone (IBZ) is used. To estimate the lattice thermal conductivity of PNRs, we 
use the phonon Boltzmann transport equation as implemented in the so-called 
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Figure 7 | The calculated electronic transport coefficients at 300 K as a function of chemical potential for the H-passivated 7, 8, 9-ZPNR: (a) the 
electrical conductivity, (b) the power factor, (c) the electronic thermal conductivity, and (d) the ZT value. 
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ShengBTE code'^^"*''. For the calculations of the second-order and third-order 
interatomic force constant matrix, we use a 1 X 1 X 6 and 1X1X3 supercell, 
respectively. The interactions up to the third nearest neighbors are considered when 
dealing with the anharmonic one. The g-point grid and scale parameter for Gaussian 
smearing are set as 1 X 1 X 35 and 1.0, respectively. 
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